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First principles calculations are used to investigate electronic band structure and vibrational 
spectra of pnictogen substituted ternary skutterudites. We compare the results with the prototypical 
binary composition CoSbs to identify the effects of substitutions on the Sb site, and evaluate the 
potential of ternary skutterudites for thermoelectric applications. Electronic transport coefficients 
are computed within the Boltzmann transport formalism assuming a constant relaxation time, 
using a new methodology based on maximally localized Wannier function interpolation. Our results 
point to a large sensitivity of the electronic transport coefficients to carrier concentration and to 
scattering mechanisms associated with the enhanced polarity. The ionic character of the bonds is 
used to explain the detrimental effect on the thermoelectric properties. 



PACS numbers: 72.20.-i,63.20.D-,71.15.Mb 

I. INTRODUCTION 

Thermoelectric materials with the filled skutterudite 
structure are considered to be a prototypical realization 
of the phonon-glass electron-crystal paradigm (PGEC) 
proposed by Slack.^ Indeed many compositions in this 
structural family exhibit low thermal conductivities (k 
0.5 - 5 Wm-^K-i), Seebeck coefficients (S) from -200 
jiV/K to 200 jiV/K., and electrical resistivities (p) in 
the range from 10~^ to 10~^ • cm at room tempera- 
ture, depending on doping levels.^ Their figure of merit 
ZT (ZT = TS^ I pk is used to characterize the material's 
performance^ reaches values in excess of 1.4 at high tem- 
perature in the bulk form.^ Skutterudites have been in- 
vestigated for use in high-reliability thermoelectric mod- 
ules designed for space applications,^ owing to their good 
thermal stability and mechanical strength throughout 
the operating temperature range. Mechanical strength 
is of particular importance in automotive and household 
applications,'^ where stress due to repeated thermal cy- 
cling is a major engineering challenge. The chemical ro- 
bustness and stability of the skutterudite crystal struc- 
ture allows for composition modifications across a wide 
chemical space, which in turn provides freedom in opti- 
mizing electronic and thermal transport properties. In 
this work, we explore one such variation: heterogeneous 
pnictogen substitution. 

The conventional cubic unit cell of a binary skutteru- 
dite such as CoSba (four formula units, space group n. 
204) consists of a simple cubic transition metal (M=Co) 
sublattice intertwined with square rings formed by co- 
valently bonded pnictogen (X4=Sb4) ions and oriented 
along (100), (010) and (001) crystallographic directions. 
Each transition metal sits at the center of a distorted 




FIG. 1. (Color onhne) Rhombohedral (R3) unit ceh of PSTS 
such as C0X1.5Y1.5 compound containing C0X3Y3 octahedra. 
Co centered octahedra linked by nearly rectangular X-Y 4- 
member rings (light gray lines), a characteristic feature of 
ternaries skutterudites. Black lines represent the unit cell 
(figure produced with CrystalMaker) . 



pnictogen octahedron. In general, the six pnictogens 
share nine electrons with the transition metals and two 
other electrons with the two nearest pnictogen ions. 
Charge balance constrains the transition metal atom to 
have 9 electrons {d^s'^) , thus leaving a limited choice of 
binaries with Co, Rh or Ir as M and P, As, Sb as X, in 
the absence of a filler ion. Two large voids, that remain 
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empty in the binary skutterudite cell, can be occupied 
by large ions (fillers).^ Substitutions and filling have a 
strong effect not only on the lattice thermal conductiv- 
ity but also on the electronic band structure and asso- 
ciated transport properties. This was pointed out both 
experimentallyP and by first principles band structure 
calculations.^™^ 

From a thermoelectricity point of view, binary skut- 
terudites have comparatively large thermal conductivity, 
k. However, filling dramatically lowers the lattice compo- 
nent of /c, (kL) and enhances ZT.^ 22 31 ^Hoying on the 
transition metal site has also been explored as a strategy 
to decrea se therma l conductivity and control electronic 
transport .f^^^MEll However, the effect of chemical substi- 
tution on the pnictogen site remains largely unexplored. 

In this work we focus primarily on recently synthe- 
sized CoGei.sSi.s, GeSni.sTei.s and CoGei.sTei.s where 
the substitution is occurring on the pnictogen site of the 
prototypical CoSba. We call these materials pnictogen 
substituted ternary skutterudites (PSTS). In order to ob- 
tain a complete comparison we also studied CoGei.sSei.s, 
CoSni.sSi.s, and CoSni.sSei.s. PSTS are experimentally 
observed to have a significantly lower thermal conductiv- 
ity than CoSbs ,1^^^ and thus are attractive to be inves- 
tigated as potential thermoelectric materials. Until now 
the features of the band structure and of the phonon dis- 
persion have not been investigated in detail. 

This paper is organized as follows. In Sec. II we 
briefly discuss the first principles methodology used to 
compute the electronic and transport properties as well 
as the phonon dispersion. Sec. Ill is devoted to the main 
results on the structural features. Sec. IV contains the 
discussion on the electronic bands , Sec. V presents the 
phonon dispersions, and Sec. VI discusses the electronic 
transport coefficients. In Sec. VII we draw our conclu- 
sions. 



II. METHODOLOGY 

All presented data are obtained by ab initio calcu- 
lations within DFT formalisnP^'^ using the Per dew- 
Zunger LDA exchange-correlation energy 

functional.ElEll 

The effect of the core electrons is treated within the pseu- 
dopotential approach with both ultrasoft (Co, S),'Sand 
separable norm-conserving (Ge, Sn, Te) pseudopoten- 
tials. Plane- wave basis was employed for the expansion 
of the valence electron wave functions and charge den- 
sities with the kinetic-energy cutoffs of 30 Ry and 240 
Ry, respectively. All calculations were performed using a 
4x4x4 Monkhorst-Pack k-point mesh to sample the Bril- 
louin zone. All internal atomic coordinates were relaxed 
within Broyden-Fletcher-Goldfarb-Shanno method until 
the forces on the nuclei were below 10~^ a.u. The theoret- 
ically optimized lattice provides a residual stress smaller 
than 5.8 KBar. 

Phonons were computed using density functional per- 
turbation theory (DFPT).^^ The dynamical matrix was 



Fourier interpolated on a fine q-point mesh starting from 
a 2x2x2 grid. All calculation were performed with the 
Quantum-ESPRESSO software.l^SI 

Electronic transport coefficients are computed with- 
ing the BOLTZWAN code^^ using the Boltzmann trans- 
port equation (BTE) in the constant scattering time ap- 
proximation. Our methodology differs from other ap- 
proaches (see for instance Ref. |48| in that we employ 
maximally localized Wannier functions (MLWF, Ref. |49|) 
to map the first principles electronic structure on a tight- 
binding model and obtain band derivatives following fol- 
low the work of Yates et al. (Ref. 1501). The method is 
not sensitive to band crossings and provides an efficient 
way to integrate Fermi velocities over the Brillouin zone. 
The computation of MLWFs has been performed within 
the WANNIER90 package^^ using the Bloch states ob- 
tained with the Quantum-ESPRESSO distribution. Rel- 
evant procedures for obtaining the band derivatives are 
described in the Appendix and examples of MLWF for 
CoSbs are shown in Fig. ^ We used also the BOLTZ- 
TRAP package^^ of Madsen and Singh to compare with 
previous calculations. 



III. STRUCTURAL FEATURES 

The two main structural units in prototypical CoSba 
are transition metal centered pnictogen octahedra and 
pnictogen rings. In PSTS the symmetry decreases with 
respect to CoSba and two different kinds of octahedra and 
rings can be identified. The structure of the pnictogen 
rings is known to have a strong influence on electronic 
bands, phonons, and consequentl y transpor t properties 
of binary and filled skutterudites j -^^^ ' ^-^ ' ^ ^l^ The typical 
PSTS structure, MX1.5Y1.5 (space group n. 148) is de- 
rived from the binary counterpart by a substitution of 
the pnictogen ion with a pair of elements from 14 (Ge, 
Sn) and 16 (S, Se, Te) groups. The stoichiometry is pre- 
served but heterogeneity is introduced in the rectangular 
rings in which the two different ions are opposite (trans) 
to each other. The rhombohedral primitive cell contains 
32 atoms and can be described as a corner sharing oc- 
tahedral network that contains two non-equivalent Co- 
sites (2c and 6f Wyckoff positions, 2c along the diago- 
nal of the cube), two non-equivalent X-sites (6f and 6f), 
and two non-equivalent Y-sites (6f and 6f). Each tilted 
octahedron is formed by group 14 and 16 ions ordered 
in alternating layer perpendicular to the [111] direction. 
In general, the pattern of Co off-center displacements is 
such that the structure is centrosymmetric. For all the 
three compounds we have analyzed, Co(2c) is off center in 
their respective octahedra and displaced along the [111] 
direction toward the smaller surrounding ions; Co(6f) is 
also slightly displaced toward the smaller ions (this true 
for all cases except CoGei.sSei.s where the covalent radii 
are very similar) but in more complex pattern compatible 
with the symmetry. The octahedral units are deformed 
and tilted (a+a+a+ in Glazer notation). The tilting is 
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TABLE L Cystal structure of PSTS C0X1.5Y1.5 after relaxation of all internal degrees of freedom. The symmetry is R3 (space 
group 148); experimental data (from Refs. [29l l34l [35l I37H39)) are between parenthesis. CoSbs has higher symmetry (space 
group 204) but is treated as R3 for easier comparison. 





CoGe 


1.5S1.5 


CoGei.sSei.s 


CoGei 


.5Tei.5 


CoSni.sSi.s 


CoSrii.sSei.s 


CoSni 


.5Tei.5 


CoSba 


ai. (A) 


7.888 


(8.010) 


8.186 


8.622 


(8.699) 


8.311 


8.610 


9.023 


(9.122) 


8.972 (9.038) 


a (deg) 


89.90 


(89.94) 


89.83 


89.95 


(89.99) 


89.87 


89.98 


89.97 


(90.06) 


90.0 


Co (2c) X 


0.258 


(0.258) 


0.251 


0.243 


(0.247) 


0.267 


0.260 


0.253 


(0.250) 


0.25 


Co (6f) X 


0.258 


(0.262) 


0.253 


0.249 


(0.249) 


0.262 


0.260 


0.255 


(0.250) 


0.25 


Co (6f) y 


0.762 


(0.755) 


0.753 


0.745 


(0.745) 


0.773 


0.764 


0.756 


(0.750) 


0.75 


Co (6f) z 


0.754 


(0.750) 


0.752 


0.747 


(0.749) 


0.758 


0.755 


0.751 


(0.750) 


0.75 


Xa (6f) X 


0.999 


(0.000) 


0.998 


0.996 


(0.995) 


0.001 


0.999 


0.998 


(0.998) 


0.000 


Xa (6f) y 


0.335 


(0.336) 


0.327 


0.318 


(0.318) 


0.333 


0.328 


0.321 


(0.319) 


0.334 (0.335) 


Xa (6f) z 


0.151 


(0.148) 


0.158 


0.167 


(0.166) 


0.149 


0.156 


0.165 


(0.162) 


0.159 (0.158) 


Xb (6f) X 


0.499 


(0.498) 


0.500 


0.501 


(0.501) 


0.499 


0.500 


0.501 


(0.500) 


0.5 


Xb (6f) y 


0.835 


(0.836) 


0.827 


0.818 


(0.829) 


0.834 


0.828 


0.821 


(0.823) 


0.834 (0.835) 


Xb (6f) z 


0.349 


(0.350) 


0.341 


0.332 


(0.338) 


0.351 


0.343 


0.335 


(0.337) 


0.341 (0.342) 


Ya (6f) X 


0.000 


(0.001) 


0.00 


0.999 


(0.001) 


0.001 


0.001 


0.000 


(0.001) 


0.00 


Ya (6f) y 


0.344 


(0.347) 


0.344 


0.345 


(0.346) 


0.337 


0.328 


0.339 


(0.338) 


0.334 (0.335) 


Ya (6f) z 


0.849 


(0.856) 


0.850 


0.851 


(0.854) 


0.840 


0.843 


0.845 


(0.845) 


0.841 (0.842) 


Yb (6f) X 


0.502 


(0.505) 


0.503 


0.505 


(0.501) 


0.501 


0.502 


0.503 


(0.503) 


0.5 


Yb (6f) y 


0.844 


(0.846) 


0.844 


0.845 


(0.842) 


0.837 


0.838 


0.839 


(0.841) 


0.834 (0.835) 


Yb (6f) z 


0.650 


(0.646) 


0.649 


0.648 


(0.652) 


0.659 


0.657 


0.655 


(0.655) 


0.659 (0.658) 



established to form the bonds of the two non-equivalent 
four- member rings involving Y a and or with and 
in the PSTS structure and involves a doubling of the 
unit cell with respect an ideal ReOa network. Shorter 
bonds are formed along a preferred cartesian direction 
and, to accommodate the rigidity of the octahedra, longer 
bonds result in one of the perpendicular directions. The 
relative length of these bonds determines the deviations 
from the ideal square shape (Oftedal's law) of the pnic- 
togen rings. In PSTS such deviations are larger than 
in CoSba since the bonds have additional ionicity that 
tends to decrease the interatomic distances (Schomaker- 
Stevenson rule). The dihedral angle in the rings changes 
from 90.0° in CoSba to smaller values ranging from 81.7° 
to 89.8° degrees for all the compounds except CoSni.sSi.s 
and CoSni.sSei.s. Our computed structural parameters 
are given m Tab. l and are within 2% of the experimental 
data.^^l^^l^^l^^Mrxhe lattice parameter correlates well 
with the covalent radii of the main group elements and 
the cell remains pseudo-cubic with small rhombohedral 
angles. Our data shows the expected correlation between 
the lattice parameter and the size of the substitution 
atoms on both pnictogen sites. For example, among the 
Ge-substituted compounds CoGei.sSi.s has the smallest 
lattice size while CoGei.sTei.s has the largest and the 
same trend also appears in the other substitution site. 



IV. ELECTRONIC STRUCTURE AND 
TRANSPORT 

We calculate the electronic band structures (Fig. [2| 
and compare them with the one of CoSba (under equiv- 
alent symmetry representations) in order to investigate 



the effect of the pnictogen substitution. In all cases the 
valence bands consist of three separate manifolds. The 
lowest two are primarily derived from the unmixed s 
states of two pnictogen types, and the splitting of the 
bands observed in PSTS is due to the different chemical 
nature and electronegativity of the pnictogen ions form- 
ing the rings. By comparison, in CoSbs the Sb-5 states 
contribute one single manifold. Both top valence and 
bottom conduction bands consist primarily of a mixture 
of Co d states and pnictogen p states with the majority 
of d states lying below the top of the valence band. 

Although the value of the computed band gap depends 
on the type of the exchange-correlation functional used, 
| 13 | 20|54 | our cases the direct gap is 2-3 times larger 
than in CoSbs (it ranges from 0.41 eV in CoSni.sSei.s to 
0.61 in CoGei.sSi.s). For comparison our calculations for 
CoSbs give an energy gap of 0.22 eV in our DFT LDA 
calculations, while the experimentally measured values 
exhibit a wide variation, i ^^ l ^^tt eij 

Several effects contribute to the change in the band 
gap, mainly the tj^-e* derived manifold splitting and the 
flattening of the band dispersion in PSTS induced by 
the more ionic bonding. Skutterudite systems typically 
possess a single band that disperses away from the tj^ va- 
lence manifold and reaches its maximum at F point (the 
highest occupied band). This band controls the lower 
edge of the energy gap and is important due to its role 
in transport in p-type materials because it provides car- 
riers with small effective mass. In CoSni.sTei.s the top 
of the valence band is about 170 meV above the low ly- 
ing d-bands. This separation increases in CoSni.sTei.s 
(220 meV), in CoGei.sSei.s (250 meV), and in the other 
PSTS, reaching values higher than in CoSbs (370 meV). 
The second higher energy valence band (from the 
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FIG. 2. (Color online) First principles band structure of C0X1.5Y1.5 with (X,Y) = (Ge,S), (Ge,Se), (Ge,Te), (Sn,S), (Sn,Se), 
(Sn,Te) (dashed) compared to binary CoSba. CoSba has been represented with R3 symmetry for comparative purposes. 



manifold) of PSTS have a multivalley character with 
heavy effective masses. Particularly in CoGei.sTei.s and 
in CoSni.sTei.s the top of the valence band is relatively 
close in energy to the bands below it; if this energy dif- 
ference could be further reduced the contribution from 
heavier carriers would enhance the p-type Seebeck coeffi- 



cient favoring the thermoelectric performance. For com- 
parison, in La filled CoSba the first heavy valence band 
is about 70 meV below the top of the valence band due 
to an interaction between filler /-states and the highest 
valence band.^^ 

In order to investigate the effects of ternary substitu- 
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tion on transport properties, we first evaluate the inverse 
of the hole effective mass tensor in the Wannier represen- 
tation (see Appendix). The inverse of the effective mass 
is then defined as an average of the diagonal elements of 
the tensor, 1/m* = ^Xl^l/^i- The corresponding val- 
ues are 0.196me, 0.169me and 0.134me for CoGei.sSi.s, 
CoGei.sTei.s and CoSni.sTei.s respectively, where mg 
is the electron mass. The se value s are larger than re- 
ported ^Omrrie for CoSbs .I^^I^^ In p-type PSTS sam- 
ples higher effective masses of carriers are presumably 
responsible for the larger Seebeck coefficient values ob- 
served experimentally. We find that the dispersion of the 
top valence band is also affected by the pnictogen sub- 
stitutions; it is more parabolic than in CoSba although 
it also exhibits linear character of the dispersion close to 
r. This linearity has been suggested in earlier work to 
affect hole transport and deviate from traditional semi- 
conducting behavior. 

The lowest conduction energy levels also exhibit new 
features in PSTS. Several non-equivalent minima in F — L 
and T — X directions can provide pockets of carriers with 
large effective masses upon n-type doping. This effect is 
also due to the decreased dispersion of pnictogen p bands 
due to stronger ionicity. 

We derive the electrical conductivity and the Seebeck 
coefficient by solving the Boltzmann transport equation 
(BTE) in the constant relaxation time (r) approxima- 
tion. We assume r = 10 fs in this work, which is com- 
monly used for studying semiconductors.^^ This is an ar- 
bitrary choice since the scattering time for the PSTS is 
not known but it allows to establish the trends associated 
with band structure effects. Within the constant relax- 
ation time approximation the Seebeck coefficient does not 
depend on r and computed values can be compared di- 
rectly with experimental data as we have done in Fig. |6j 
[7l and|8l 
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FIG. 3. (Color online) Seebeck coefficients of ternary skut- 
terudites at 300K as a function of the electron chemical po- 
tential (Fermi energy). 

At the room temperature the Seebeck coefficient of 
CoGei.sSi.s ranges from 39 to 258 \iV jK for p-type dop- 
ing in the range of 10^^ to 10^^ electron/cm^. Since 
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FIG. 4. (Color online) Electrical conductivity of ternary skut- 
terudites at 300K as a function of the electron chemical po- 
tential (Fermi energy). 
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FIG. 5. (Color online) Power factor of ternary skutterudites 
at 300K as a function of the electron chemical potential (Fermi 
energy). 



the experimental value of the carrier concentration is 
not available for the samples under investigation we per- 
formed computations of the Seebeck coefficient in a range 
of possible carrier concentration by varying the position 
of the electron chemical potential. For carrier concentra- 
tion of 10^^ holes per cm^ our results agree with exper- 
imentally reported thermopower of CoGei.sSi.s. In the 
case of CoGei.sTei.s and CoSni.sTei.s, following the ex- 
perimental findings, we tested n-type doping in the limits 
of 10^^ to 10^° electron per cm'^. At room temperature S 
varies between -646 to -257 fiV/K for CoGei.sTei.s and 
from -695 to -307 iiV/K for CoSni.sTei.g. While com- 
puted values for CoSni.sTei.s are in reasonable agree- 
ment with experimental data at a carrier concentration 
of 10^^ per cm^, our results differ from the experimental 
data for CoGei.sTei.s at low temperatures. The See- 
beck value reaches its minima of nearly -800 jiV/K at 
115 K, exhibiting an apparent dip and a subsequent in- 
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crease decrease in the magnitude. It may be tempting 
to explain such a trend reversal by the bipolar effect, i.e. 
decrease of S due to thermal activation of minority car- 
riers across the band gap. However, we argue that this 
feature derives only from the electronic structure of the 
valence manifold at the experimental carrier concentra- 
tion. The first reason is that the bipolar effect typically 
sets in at temperatures where ksT is comparable to the 
band gap, and in PSTS the band gap is comparably large. 
The experimental values of S for CoGei.sSi.s, whose 
atomic and electronic structure is similar, exhibit the 
expected trend of increasing S with temperature. The 
second and most compelling reason is the appearance of 
such a non-monotonic dip feature in the computed See- 
beck coefficient temperature dependence of CoSni.sTei.s 
(Fig. [s]) at n=10^^ cm~^. From the computed electronic 
band structure of CoSui.sTei.s (Fig. [2| we readily con- 
clude that this non-monotonic feature is due to the inter- 
twining and non-dispersive character of the valence band 
manifold. The deviation from the experimental behav- 
ior at higher temperatures is likely due to the the varia- 
tion of the actual carrier concentration with temperature 
and the inaccuracy of the constant r approximation.^^ 
At low temperatures impurity states may also account 
for the observed strong dependence of transport prop- 
erties on temperature.^^ For CoGei.sSei.s, CoSni.sSei.s, 
and CoSni.sTei.s we find values of S to be between 200 
and 400 fiV/K for n-type doping and between 400 and 
600 iiV /K at room temperature. The general trend, as 
shown in Fig. [3] is that the Seebeck coefficient in all six 
PSTS increases substantially with respect to CoSbs both 
for p-type and n-type doping; this reflects the decreased 
band dispersion. 

Experimental values of electronic resistivities at the 
room temperature are reported 30.6 l^cm,^^ for p-type 
CoGei.sSi.s , 5.1 (^cm,^^, for n-type CoGei.sTei.s, and 
0.33 Qcm,^for n-type CoSni.sTei.s. The theoretical re- 
sults for conductivity at room temperature span over sev- 
eral orders of magnitude depending on the carrier con- 
centration. Our approach is to determine the carrier con- 
centration as the one that produces the best match the 
temperature dependance of the thermopower to experi- 
mental measurements. 

Selecting this doping level, we compute the room tem- 
perature electrical conductivity that is about two orders 
of magnitude larger than experimental values. It must be 
noted that our methodology reproduces the experimental 
results in wide range of temperature for the conductivity 
of CoSba assuming r = 2.5 10""*^^ (Ref. "2T) when we use 
the experimentally determined carrier concentration val- 
ues. Within the constant relaxation time approximation, 
where thermopower S is independent of r, this approach 
provides a way to separate possible contributions to the 
discrepancy between theory and measurements. We can 
reasonably conclude that the features of the electronic 
structure alone are only partially responsible for much 
larger electrical conductivity with respect to experiment. 
Three quantities contribute to the electronic conductiv- 



ity: effective masses, carrier concentration and scatter- 
ing time r. The reasons for the discrepancy between 
experiment and theory could include inaccurate carrier 
concentrations and/or an anomalously short carrier life- 
time. Impurity phases and defects in the experimental 
samples may also contribute to the discrepancy with the 
computed results. 

In order to evaluate the potential of PSTS as active 
materials in thermoelectric devices, we compare their 
performance with the well-known CoSba material. All 
PSTS have lower electronic conductivity than in a wide 
range of doping levels, as shown in Fig. [4j Since the value 
of r is taken to be the same, this reflects the larger band 
gap, decreased band dispersion and larger carrier effec- 
tive masses. In Fig. [5] we show the full power factors of 
all compositions as a function of doping level, and these 
results show a noticeably lower power factor for PSTS as 
compared with CoSbs in the p-type region and most of 
the n-type region. We conclude that, in the electronic 
transport aspect, PSTS are not likely to surpass the per- 
formance of CoSbs-based systems, particularly for p-type 
materials, assuming the same carrier lifetimes. Further- 
more, as we discuss below, carrier lifetimes in PSTS are 
likely reduced by the enhanced ionicity. 

In an ideal crystal the scattering time includes con- 
tributions from the electron-phonon coupling, the larger 
contribution associated with deformation potential and 
Frohlich scattering. The enhanced ionicity in PSTS sug- 
gest to consider effects associated to the Frohlich inter- 
action. We have qualitatively analyzed this contribution 
by evaluating the mode-resolved Born effective charges, 
defined by 

^a(^.k) = = (1) 

to estimate for the polarization arising from the vibra- 
tional displacements and, consequently, the strength of 
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FIG. 6. (Color online) Temperature and doping dependence 
of the Seebeck coefficient of CoGei.sSi.s compared to experi- 
mental intrinsic data."^ 
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FIG. 7. (Color online) Temperature and doping dependence 
of the Seebeck coefficient of CoGei.sTei.s compared to exper- 
imental intrinsic data^^. 
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FIG. 8. (Color online) Thermopower of 

CoSni.sTei.sfbottom) compared to the experimental in- 
trinsic data.^^ Note the dip at low temperature as discussed 
in the text. 



electronic scattering."^ 

Due to the smaller primitive cell and weak ionicity, 
CoSba has a few (7) vibrational frequencies that exhibit 
non-zero Za,{y^^ k) at the Brillouin zone center (see Sec. [V| 
for the phonon dispersions). In the low frequency region 
below 120 cm~^ the mode resolved effective charges are 
less than 1 and the vibrational modes do not effectively 
scatter electrons. More significant scattering is expected 
when modes above 250 cm~^, with ZcK(a;, F) 2^ 8, become 
active. In PSTS the situation is quite different: we com- 
puted, in fact, many "polar modes" with 2-3 time larger 
Za, than in CoSbs. These modes are distributed across 
the entire frequency spectrum. This indicates that the en- 



hanced polar scattering, especially at low frequency, may 
affect strongly the electrical conductivity, as compared to 
CoSba. It is important to notice that the polar scattering 
contribution affects the thermal conductivity as well. 



V. PHONONS 

Full first principles phonon dispersion for filled an d un- 
filled skutterudites were studied by Feldman et al.j "'^^ * "'^^ * 
Ghosez et al.,1^ and Wee et al.l^ The vibrational spec- 
trum of PSTS is an essential starting point to understand 
the role of the chemical substitutions in PSTS and de- 
velop models for the low thermal conductivity observed 
in these materials. We present here the vibrational dis- 
persions at the theoretically optimized structural param- 
eters. For comparison, in CoSba there are two main man- 
ifolds associated, respectively, with the vibration of the 
transition metal (between 250 and 300 cm~^) and of the 
pnictogens (below about 200 cm~^). In CoGei.sSi.s (Fig. 
|9| the comparable masses of Co and Ge result in the for- 
mation of vibrational modes that are mixed in character. 
The dispersion shows an additional manifold associated 
mainly with sulphur vibration above 350 cm~^. The mo- 
tion of Co contributes across all frequencies with a larger 
contribution near 300 cm~^. The frequency of the low- 
est optical mode (mainly Ge) at F is at about 100 cm~^ 
only slightly higher then the Sb-modes in CoSbs. Simi- 
lar features are observed in the dispersion of CoSni.sSi.s 
(not shown) where, of course, the Sn-derived modes ex- 
tend to lower frequencies (about 75 cm~^). The phonon 
dispersions of CoGei.sTei.s (Fig. 10) and CoSni.sTei.s 
(Fig. 11) exhibit two manifolds, similar to CoSba: the 



highest manifold is mostly from Co motion. The lowest 
frequency optical modes are at 65 cm~^ in CoGei.sTei.s 
and 50 cm~^ in CoSni.sTei.s. This is the frequency re- 
gion where modes from filler atom vibrations are found 
in BaCo4Sbi2 (Ref. 21) and may point to the phonon 
scattering channel responsible for the low thermal con- 
ductivity. 

The group velocity of acoustic modes near F deter- 
mines the thermal conductivity and, in our calculations, 
correlates with the mass of the specific pnictogen sub- 
stituted ions. It is interesting to notice that the sound 
velocities in CoSni.sTei.s are very similar to those of 
CoSba. In other PSTS we found values higher than those 
of CoSba. Based on these results and the overall phonon 
dispersions, it is reasonable to argue that scattering phe- 
nomena differ substantially between CoSba and PSTS 
probably due to the different character of the bonding 
in the rings. Phonon dispersions alone cannot explain 
the low thermal conductivity values of observed experi- 
mentally for PSTS. More work in the direction of under- 
standing the anharmonic scattering in these materials is 
required. 
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TABLE 11. Transverse effective charges, Z*, computed with density functional perturbation theory. The full tensors are used 
to compute the electron-phonon polar scattering contribution but only |TrZ* is reported here. 



X = Ge,Y = S X = Ge,Y = Te X = Sn,Y = Te X = Y = Sb 



Co (2c) -5.140 -6.021 -5.914 -6.678 

Co (6f) -4.895 -5.900 -5.810 

Xa +3.227 +3.463 +3.506 +2.229 

X B +3.195 +3.462 +3.480 

Ya -0.035 +0.439 +0.334 

Yb -0.006 +0.430 +0.353 



VI. CONCLUSIONS 

We discussed structural aspects, electronic structure 
and transport, and phonon dispersions of pnictogen sub- 
stituted ternary skutterudites (PSTS). These materials 
are potentially interesting for thermoelectric applications 
due to the exhibited low lattice thermal conductivity. 
Unfortunately the electronic transport is not as favorable 
because of the low electrical conductivities. 

We justified the large Seebeck coefficients by analyzing 
the electronic bands structures: a decreased dispersion 
compared with CoSba as well as a multivalley character 
with heavy carrier effective masses. The values of elec- 
tronic conductivity are lower than for CoSba and have 
a strong dependence upon carrier concentration. We ex- 
plored the upper limits on the power factor of PSTS in 
a wide range of carrier concentrations and found that 
they are unlikely to surpass those of CoSba. More ef- 
fort should be invested in understanding the reasons for 
low measured values and find way to increase electronic 
conductivity in these materials. 
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VIII. APPENDIX: BOLTZMANN TRANSPORT 
FROM WANNIER FUNCTIONS 
INTERPOLATION 

Prediction of electronic transport properties, using the 
Boltzmann transport equation (BTE), depends on the 
ability to accurately compute and integrate band deriva- 
tives over the Brillouin zone. Usually this is achieved 
by fitting the electronic band to a smooth curve and 
performing numerical derivatives, an approach that is 
sensitive to band crossings. The Wannier representation 
of the electronic structure^^ 66 68 provides an optimized 



tight-binding model whose Hamiltonian can be directly 
differentiated to compute band velocities and effective 
masses. ^^-^ An additional advantage of the approach is 
the possibility to separate the role of individual bands 
or band manifolds by projecting on minimal subspaces 
containing the most relevant degrees of freedom, using 
the disentanglement procedure.!^ For transport proper- 
ties only a certain subset of the Bloch states near the 
Fermi level is relevant. In this work we used Maximally 
Localized Wannier Functions (MLWF) to derive the nec- 
essary ingredients for the BTE in the constant scattering 
time approximation. As a side product we obtained a de- 
scription of the bonding states in terms of MLWF (Fig. 
12). 



Prototypical CoSba is a semiconductor with two iso- 
lated valence manifold of 12 and 36 bands respectively 
and a conduction manifold that consists of infinite num- 
ber of entangled states above a LDA energy band gap of 
the order of ~ 0.22 eV. The lowest manifold of 12 valence 
states is mainly formed by Sb s-states. The top 36-band 
manifold is constructed iteratively from the initial guess 
of the atomic Sb s- and p- and Co d-states. Starting 
with a combination of on-site Co d-states and Sb s- or p- 
states we have converted original combinations of atomic 
orbitals to a well localized set of Wannier functions with 
spreads in the range of 1.5-6.65 A^. Among all 48 valence 
states one can distinguish 12 Co states of t2g symmetry, 
12 Sb-Sb bonding states and 24 Co-Sb bonding states 
(Fig. 12). To construct Wannier states for the conduc- 
tion manifold for CoSba we choose Bloch states in the en- 
ergy range of 3.2 eV above the Fermi level. MLWF states 
were obtained by iterative convergence starting with the 
initial guess of 24 gaussian-type orbital states placed 1/4 
of the Co-Sb bond length away from Co atoms along each 
of the 24 Co-Sb bonds. Using a similar approach we have 
also determined the basis of MLWFs for the PSTS sys- 
tems. 

Given the basis of MLWFs we can express matrix ele- 
ments of the Hamiltonian in terms of Wannier functions. 
Matrix elements of a periodic operator O between Wan- 
nier states n and m are written as 0^^{R) = {nO\0\mR) 
The matrix element of the Hamiltonian at an arbitrary k 
point in the k-space (Ref 42.) can be obtained by inverse 
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FIG. 9. (Color online) Calculated phonon dispersion and atom projected vibrational densities of states of CoGei.sSi.s- 



Fourier Transformation (FT) interpolation 
= ^e'''^{nO\H\mR) 

R 

. Due to the strong localization of MLWFs the Hamilto- 
nian in the Wannier basis is sparse and one does not need 
the original k-point mesh (used to construct the Wannier 
states) to be dense to obtain convergence for an arbi- 
trary k point. For large systems, FFT scales much faster 
{0{N\og{N))) compared to the scaling of the eigenvalue 
problem ( 0{N^)). The right hand side of the last equa- 
tion can be differentiated analytically with respect to k 
to obtain the matrix elements of the velocity operator. 



^ R 

As a last step, a rotation to the original set of the Bloch 
states is performed. This however, requires matrix mul- 
tiplication of only very small matrices of M x M size, 
where M is the number of Wannier states. As a result, 
this interpolation scheme is faster than the direct solu- 
tion of the eigenvalue problem, but with the additional 
complexity of the initial Wannierization. It also resolves 
a number of difficulties associated with band crossings 
and avoided crossings which persist in traditional inter- 
polation schemes. 



10 




FIG. 10. (Color online) Calculated phonon dispersion and atom projected vibrational densities of states of CoGei.sTei.s. 
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